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Abstract 

g ■ We investigated the effect of on-site Coulomb interactions on the structural and magnetic ground state 

, of the fullerene C20 based on density-functional-theory calculations within the local density approximation 

Q ' plus on-site Coulomb corrections (LDA-i-f7). The total energies of the high symmetry (Ih) and distorted 

(D^fi) structures of C20 were calculated for different spin configurations. The ground state configurations 
were found to depend on the forms of exchange-correlation potentials and the on-site Coulomb interaction 
ly^ I parameter U, reflecting the subtle nature of the competition between Jahn-Teller distortion and magnetic 

CN ' instability in fullerene C2o- While the non-magnetic state of the distorted D^^ structure is robust for small 

^ ■ U, a magnetic ground state of the undistorted 1^ structure emerges for U lai^ger than 4 eV when the LDA 

Q , exchange-correlation potential is employed. 
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I. INTRODUCTION 



Since the discovery of fuUerene Ceo UJ] various research activities on alternative fuUerene struc- 
tures of carbon materials have flourished due to their noble physical properties and potential appli- 
cations as next generation electronic devices. The smallest fuUerene C20 has been actively studied 
after its first successful production by Prinzbach et al. ||2|]. Although there have been many the- 
oretical predictions on its intriguing properties such as superconductivity, [sj] vibrational modes, 
l^lsj] and transport 1^7,8] along the line of fuUerene-based molecular electronics, [S,!!^ C20 has 
been a challenge for theorists since long before its experimental synthesis. For instance, the stable 
structural configuration of C20 has not been clearly understood. The effort to determine the stable 
structure and relative stability of C20 isomers, e.g., ring, bowl, and cage, has raised a question 



IJ. One 



on the contribution of electron correlation energy to the structural stability] ill [l^ 
of the difficulties in predicting the most stable C20 isomer lies in the treatment of the exchange- 
correlation (XC) energy. It was shown that density-functional-theory (DFT) calculations gave 
qualitatively different results depending on the XC functional such as local density approximation 
(LDA), generalized gradient approximation (GGA), and Hatree-Fock (HF). Even the quantum 
Monte Carlo (QMC) results were different from LDA, GGA, HF and hybrid DFT[[l3t]. 

The prediction of the ground state atomic structure of fuUerene C20 has suffered the similar 
problem. The high symmetry Ih structure of the C20 cage is known to be unstable with respect to 
the Jahn-Teller (JT) distortion. According to Galli et al. lll4ll . both LDA and GGA predicted D^^ as 
the symmetry of the most stable Jahn-Teller distorted cage structure, which is consistent with the 
tight-binding molecular-dynamics (TBMD) result by Yamamoto et al. [|8D. However, the hybrid 
DFT calculation by Saito and Miyamoto M determined to be the most stable one, while C2 or 
D^d were chosen by the HF calculations 11 1411 . Sometimes even the same kind of calculations made 
by different groups have given different predictions (see, for example, Ref. [14]). As Grossman 
et al. pointed out, this small molecular carbon system, C20, demonstrates the importance of the 
highly accurate treatment of electron correlations list]. 



Recently Lin et al. jlS . 



id ] took a different perspective to this issue. By considering the 



fuUerene C20 as a superatom with 20 valence electrons, they set up a tight binding model tak- 
ing account of the strong hybridization of ^Tr-orbitals with p^, which arises from the extreme 
curvature of the C20 cage. Expecting the important role of the on-site Coulomb interaction, they 
carried out QMC and exact diagonalization (ED) calculations for the Hubbard model hamiltonian 
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for C20. As a result, it was found that C20 with both Ih and D^d structures undergoes a magnetic- 
to-non-magnetic and non-magnetic-to-magnetic transitions for Ih and D^d, respectively, as U/t 



increases 111 5 



To investigate the effect of the on-site Coulomb interaction in connection with its ground state 



atomic and magnetic structures, we performed LDA+U and GGA+U calculations [11711 and calcu- 
lated optimized structure under the symmetry constraint of Ih and D^d- Our result showed that the 
non-magnetic D^d structure is more stable than magnetic Ih for small values of U less than 4 eV 
regardless of the choice of the XC functional, which is consistent with previous works Isi 14]. 
However, when U > A eV, the LDA+f/ calculations predict a magnetic state with /i = 2/iB to be a 
ground state, which is different from GGh+U . This result indicates that there are subtleties in the 
competition between the Jahn-Teller structural instability and the magnetic instability enhanced by 
Coulomb correlations. Further the dependence on the exchange-correlation functionals suggests 
the importance of electron correlation effects in determination of the C20 ground state. 

II. COMPUTATIONAL METHODS 
A. LDA+C/ Method 

In order to describe on-site Coulomb interactions for the carbon 2p orbitals ,which are relatively 



localized in a small-sized cluster like C20, we employed the LDA+?7 (GGA+f/) method [|l7 . 
In our LDA+f/ (GGA+f/) approach, the total energy functional is written as lllStl 



IJl. 



E-LDA+U [GGA+U) — -E'lDA(GGA) + E^, — (1) 

where Elda{gga) is the LDA (GGA) energy functional and is given by a spherically averaged 
screened-Coulomb energy U and an exchange energy J: 

= ^ E E <rnn^^^' + ^ " '^«) E <mKm' , (2) 

a amm' a a^m^m' 

where o is the spin index, ol = (Up) with the site index i, the angular momentum quantum number 
/, and the multiplicity number of the radial basis function p. is an eigenvalue of the occupation 
number matrix. Here the U- and the J-values are assumed to depend only on the index a, but to 
be independent of the azimuthal quantum number m, which is regarded as a simplification of the 
Hartree-Fock theory by using a spherical average. The double-counting term E^^ can be written 
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as 



where 



Edc =IY1 UaNaiN^ - 1) - ^ E E - 1 

Km and N^ = Nl + Ni. Therefore, Eu = El 

Eu = \ E(^" - E {Tr«) - Tr(«)} 



(3) 



E^^ reads 



(4) 



and it is the term that describes the on-site Coulomb interactions. The effective [/-value, defined 
hyU = Ua — Ja, can be treated as a numerical parameter, and it is clear from Eq. (4) that LDA+f/ 
(GGA+f/) method is equivalent to LDA (GGA) in f/ ^ limit. We used the previously proposed 



m. 



'dual' representation [|18|] . which is proven to be successful for various materials [|19l . 

We considered the on-site Coulomb interaction U for the carbon 2p orbitals, which are rel- 
atively localized in a small-sized cluster like C2o- In our LDA+f/ and GGA+f/ calculations, 
we used the previously proposed 'dual' representation [llSll . which is proven to be successful for 



various materials II 



211. 



B. Computation Details 



We performed cluster calculations based on the DFT by employing a linear-combination-of- 



23|] and Perdew- 



localized-pseudo-atomic orbitals (LCPAO) as a basis set. Ceperley-Alder [|22L 
Burke-Ernzerh [|24ll type exchange-correlation energy functional were adopted for LDA and GGA 
calculation, respectively. We used double valence orbitals as a basis set which were generated by 
a confinement potential scheme with the cutoff radius of 5.0 a.u.. TrouUier-Martins type pseudo- 
potentials with a partial core correction were used to replace the deep core potentials by norm- 
conserving soft potentials in a factorized separable form with multiple projectors. The real space 
grid techniques were used with the energy cutoff of 160 Ry in numerical integrations and the 
solution of the Poisson equation using fast Fourier transformations (FFT). In addition, the projector 
expansion method was employed to accurately calculate three-center integrals associated with the 
deep neutral atom potential with Lmax = 6 and A'rad = 4. All the DFT calculations were performed 
using our DFT OpenMX code [|25n . 
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III. RESULT AND DISCUSSION 



Figure [U shows a schematic diagram of the Hiickel's molecular orbital levels around the Fermi 
energy for the magnetic Ih and the Jahn-Teller (JT) distorted D-^d structure of fuUerene C20. In the 
highly symmetric Ih structure, the highest-occupied-molecular-orbital (HOMO) states are four- 
fold degenerate and partially filled by two valence electrons. Due to the electronic degeneracy, 
the Ih structure C20 is Jahn-Teller active and can undergo a structural distortion. Consequently 
the Jahn-Teller distorted D^d structure leads to a finite gap in between the HOMO level and the 
lowest-unoccupied-molecular-orbital (LUMO) states. Since the HOMO state is occupied by a 
singlet pair, the D^^ ground state becomes non-magnetic. On the other hand, the Ih C20 ground 
state can be magnetic due to a kind of the Hund's rule coupling present in this 'superatom', where 
the magnetic exchange energy between the HOMO electrons favors the magnetically aligned spin 
state of the undistorted Ih structure. 

The calculated density-of- states (DOS) in Fig. [2] clearly shows the different electronic structure 
of these two structures. The dotted blue and solid red lines correspond to the U = and U = 2e.Y 
calculation, respectively. The Fermi level of Ih C20 is located at the middle of a majority spin peak 
(Fig. Oa)), whereas D^d has a finite HOMO-LUMO gap (Fig. |2tb)), which is consistent with the 
Hiickel's molecular orbital pictures in Fig. \T\ While the electronic states of the JT-distorted D^d 
structure are hardly affected by the change of on-site U values from to 2 eV, the exchange split 
of the undistorted Ih structure is significantly enhanced even for the moderate on-site interaction 
of U=2 eV. Considering the enhanced exchange interactions, we tried to explore a possible con- 
tribution of the on-site Coulomb interactions to the stabilization of the magnetic Ih structure over 
that of the JT structural distortion. 

Figure |3] shows the calculated total energies by IDA+U (Fig.llta)) and GGA+U (Fig.|3];b)) as 
a function of U. For comparison, we also calculated the total energies of the non-magnetic (NM) 
phase of Ih- The total energies of FM-//j, NM-//j, and NM-D^d are represented by the dotted 
grey, dotted blue, and solid red lines, respectively. The total energy of FM-Ih is set at zero as 
a reference. For the both LDA and GGA calculations with U = 0, the NM-D^d configuration 
was found to be the ground state and the NM-J/j configuration is the highest in energy, which 
is consistent with the JT structural instability of Ih discussed above. GGA predicts the larger 
stability of the non-magnetic NM-D^d over the magnetic FM-Ih- The stabilization energy of the 
NM-Dsd over FM-4 and NM-4 is 239 and 374 meV, in LDA, whereas 821 and 998 meV in GGA, 
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, respectively. (See the Table U for detailed information.) The relative order of stability among the 
different configurations remains the same even with U = 2 eV. 

It is noted that the D-^d structure is more stabilized over 1^ in GGA calculation than in LDA 
for both U = and 2 eV; the stabilization energy is more than three times larger than in LDA. 
With a larger Coulomb interaction parameter, U = A eV, The total energy differences between 
FM-Ih and NM-D^d are dramatically reduced for both LDA+U and GGA+U results. Indeed, in 
the case of LDA+f/, the ground state is changed from the NM-D^d to the FM-//j configuration 
as shown in Fig. [31 This result implies that, when U becomes large, the magnetic energy gain 
with an enhance exchange interaction can stabilize the symmetric Ih structure against the Jahn- 
Teller distortion with a lower symmetry D^d structure. The strong on-site Coulomb interactions 
enhance the electron localization, which leads to the exchange energy gain thereby contributing to 
the stabilization of the undistorted structure. Figured shows a charge density difference between 
U = 2 and f/ = of FM-I^. The blue and red colors represent the electron surplus regions of the 
U = 2 and f/ = eV calculations, respectively. Though the magnitude of the charge difference 
is small, the difference plot of Fig. |4] demonstrates the tendency of electron localization at each 
carbon site, which arises from the electron correlation due to the on-site Coulomb interactions. The 
effect of the electron localization can affect the relatively enhanced exchange coupling among the 
degenerate HOMO states, and its magnitude seems to be comparable to that of the JT distortion. 
In Fig. [21 it is shown that the JT level spacing in Fig. [2tb) is comparable to that of the exchange 
splitting in Fig. [2ta). Since the JT distortion lifts up the orbital degeneracy of the HOMO states, 
the JT mechanism acts against the formation of magnetic moment in fuUerene C2o- This result 
is another example demonstrating the discrepancy caused by the different XC energy functionals, 
and further is consistent with the previous studies [13, [l^ which emphasized that the structural 



ground state depend on the type of XC energy functionals. 

It is interesting to compare our LDA+U and GGA+f/ calculation results with the recent 
Hubbard-model-based study by Lin et al. [[isj . Starting from an observation that the hopping 
parameter t of the C20 fuUerene should be much smaller than that of Ceo but its on-site repulsion 
U remains the same, which means that the ratio of f//t is large, they assumed that electron corre- 
lations play an important role in this system and performed ED and QMC calculation for one-band 
Hubbard model parameterized hy U/t, 

H=-t i4aCja + H.c)+UY,n^^n^l■ (5) 

<i,j>,a i 



Their results predicted that the ground state of the Ih structure changes from triplet to singlet at 
about U/t ~ 4.10, while the ground state of D^d evolves from the singlet (non-magnetic) state 
at U = to a triplet state for U/t larger than 0.5, and then transit to a singlet state again at 
U/t ^ 4.19. Those predictions based on the model analysis are in contradiction to our DFT 
calculations, where the FM state of the Ih structure is always stable relative to the NM state and 
the JT distorted D^d structure prefers a non-magnetic ground state. The main difference between 
the Hubbard model approach and our DFT calculations lies on that our DFT calculations take 
account of the exchange coupling among the electrons occupying the HOMO states. The explicit 
treatment of the degenerate HOMO states is crucial since the energetic competition between the 
exchange coupling and the JT structural distortion is crucial for the degenerate states. 

rv. CONCLUSION 

We investigated the effect of on-site Coulomb correlations on the structural and magnetic prop- 
erties of the fuUerene C20 by carrying out the LDA+f/ and GGA+U calculations. From the 
comparison of the total energies of the magnetic Ih structure and the JT-distorted D-^d structure 
for different values of the on-site Coulomb interaction U, we suggest that the exchange-coupling 
driven magnetic instability can possibly override the Jahn-Teller structural distortion at least for 
the value of U = 4 eV within the LDA+f/ calculations. Considering the subtle dependence of 
the ground state properties of C20 on the choice of XC-functional forms as well as the Hubbard 
model parameters, more elaborate calculations are required to resolve the physical picture for the 
ground state of the fuUerene C20, where the magnetic instability due to the Coulomb correlation 
effect competes with the structural distortion in the 'superatom' fuUerene C2o- 
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TABLE I: LDA+U and GGA+U total energies (in meV unit) of the Ih and D^a structures of C20 with FM 
and NM spin configurations as a function of U. FM-D^^ converges to be NM spin state in [/ = 0, 2 eV (see 
the text). The total energies are given relative to that of FM-Ij^. 



XC-type Symmetry U 


= OeV U 


= 2eV [/ 


= 4eV 


h (FM) 











LDA Ih (NM) 


135 


217 


566 


Did (NM) 


-239 


-263 


33 


h (FM) 











GGA 4 (NM) 


177 


259 


795 


Dsd (NM) 


-821 


-843 


-364 
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FIG. 1: (Color online) Hiickel's molecular orbital levels for the highly symmetric Ih structure and the Jahn- 
Teller distorted D-^d structure of €20- Up and down arrows indicate up and down spin electrons, respectively. 
The horizontal dashed line represents the Fermi level. 
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FIG. 2: (Color online) Total density-of-states (DOS) for (a) Ih and (b) D^d structure C20 fuUerene. Up and 
down pannels refer to the up-spin and down-spin states, and the dotted blue and solid red lines correspond 
to [/ = and 2 eV, respectively. The DOS plots were drawn with the Gaussian broadening of 0. 1 eV, and 
the Fermi energy is set at zero. 
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FIG. 3: (Color online) Total energy curves as a function of U. The calculated energies of NM-Ih (dotted 
blue) and NM-D^d (solid red) are plotted with respect to that of FM-//j (dotted gray) in meV unit. 



13 



FIG. 4: (Color online) Iso-surface plot of the charge density difference of FM 1^ between ?7 = eV and 
U = 2 eV. Blue- and red-colored surfaces correspond to electron surplus regions for [/ = 2 eV and [7 = 
eV, respectively, at the isovalue of 0.7 x lO^-^ e-A^'^. 



14 



